Nonequilibrium quantum criticality in bilayer itinerant ferromagnets 



O 

(N 

Oh: 

< 

(N 



CO 



So Takei 1 , William Witczak-Krempa 2 and Yong Baek Kim 2,3 
1 Max-Planck-lnstitut fur Festkorperforschung, D-70569 Stuttgart, Germany 
1 Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7, Canada 
3 School of Physics, Korea Institute for Advanced Study, Seoul 130-722, Korea 
(Dated: April 13, 2010) 

We present a theory of nonequilibrium quantum criticality in a coupled bilayer system of itinerant electron 
magnets. The model studied consists of the first layer subjected to an inplane current and open to an external 
substrate. The second layer is subject to no direct external drive, but couples to the first layer via short-ranged 
spin exchange interaction. No particle exchange is assumed between the layers. Starting from a microscopic 
fermionic model, we derive an effective action in terms of two coupled bosonic fields which are related to the 
magnetization fluctuations of the two layers. When there is no interlayer coupling, the two bosonic modes 
possess different dynamical critical exponents z with z = 2 (z = 3) for the first (second) layer. This results in 
multi-scale quantum criticality in the coupled system. It is shown that the linear coupling between the two fields 
leads to a low-energy fixed point characterized by the larger dynamical critical exponent z = 3. We compute 
the correlation length in the quantum disordered and quantum critical regimes for both the nonequilibrium 
case and in thermal equilibrium where the whole system is held at a common temperature T. We identify an 
effective temperature scale Teg with which we define a quantum-to-classical crossover that is in exact analogy 
with the thermal equilibrium case but with T replaced by T eff . However, we find that the leading correction 
to the correlation length in the quantum critical regime scales differently with respect to T and T^. We also 
note that the current in the lower layer generates a drift of the magnetization fluctuations, manifesting itself as 
a parity-breaking contribution to the effective action of the bosonic modes. In this sense, the nonequilibrium 
drive in this system plays a role which is distinct from T in the thermal equilibrium case. We also derive the 
stochastic dynamics obeyed by the critical fluctuations in the quantum critical regime and find that they do not 
fall into the previously identified dynamical universality classes. 

PACS numbers: 03.65.Yz, 05.30.-d, 71.10.-w 
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I. INTRODUCTION 

Understanding the role of nonequilibrium drives on systems 
tuned near a quantum critical point has received growing at- 
tention in recent years. One of the issues is to understand 
to what extent the concept of universality persists when sys- 
tems depart far from equilibrium. Several works considered 
universal scaling behaviour in nonlinear transport properties 
close to a superconductor-insulator quantum critical pointJr— . 
Universal scaling of current in response to an external field 
was also studied for an itinerant electron system near a mag- 
netic quantum critical point! More recently, universal scaling 
functions of the non-linear conductance and fluctuation dis- 
sipation ratios were obtained for a magnetic single-electron 
transistor near quantum criticality^. 

Complimentary works have studied the effects of cur- 
rent flow on 2D itinerant electron systems near a 
ferromagnetic-paramagnetic quantum critical point by gener- 
alizing the perturbative renormalization group to nonequilib- 
rium situations 7 -^. These works developed scaling theories 
which showed that the leading effect of the nonequilibrium 
probe is to act as an effective temperature, T e g. As in the 
corresponding equilibrium theor y 10 ' 11 , T e ff was found to be a 
relevant perturbation at the equilibrium gaussian fixed point 
and to be responsible for inducing a quantum-to-classical 
crossover similar to the one triggered by thermal fluctuations. 
A large-Teff quantum critical regime, analogous to the corre- 
sponding regime in equilibrium quantum criticality, was also 
identified. In this regime, the long-time, long-wavelength be- 
haviour of the order parameter fluctuations was found to obey 
Langevin dynamics similar to one of the dynamical universal- 



ity classes considered by Hohenberg and Halperin?. This has 
led to an important identification of a nonequilibrium dynam- 
ical universality class. 

The form of the long-time, long-wavelength dynamics 
obeyed by order parameter fluctuations depends intimately on 
the geometry of the system and on the manner in which the 
nonequilibrium perturbation is applied. For a 2D system, the 
orientation of the current flow with respect to the plane of the 
system influences the dynamics. When inversion symmetry 
is broken-, for instance with the application of current paral- 
lel to the plane, the dynamics obeyed by the order parameter 
fluctuations differs from the case where the current is applied 
normal to the plane. Introducing particle exchange by tunnel- 
coupling the system to an external reservoir also affects the 
dynamical critical phenomena. A system, initially character- 
ized by conserved order parameter fluctuations in the closed 
case, must be described by non-conserved fluctuations once 
it is opened. For itinerant ferromagnetsii this corresponds 
to a change in the dynamical critical exponent from z — 3 to 
z = 2ii The sensitivity of order parameter dynamics to system 
geometry and the form of nonequilibrium perturbation opens 
up the possibility to explore a variety of nonequilibrium dy- 
namical universality classes and contribute to the general un- 
derstanding of nonequilibrium quantum cricality. 

In this work, we consider nonequilibrium quantum critical- 
ity in a coupled bilayer system of itinerant ferromagnets tuned 
close to a spin-density wave instability. The two itinerant fer- 
romagnets are separated by a 3D insulating barrier (see FigQ]), 
and are coupled via short-ranged spin exchange interaction. 
We consider the case where there is no particle exchange be- 
tween the layers. Departure from equilibrium is achieved by 
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driving one of the layers with an in-plane electrical current. 
The resulting driftSill of the spin fluctuations, in turn, causes 
the fluctuations in the other layer to also drift, ultimately driv- 
ing the latter out of equilibrium. A steady state is established 
by tunnel-coupling the driven layer to an external bath. 

We show that the system can be characterized by two 
bosonic fields, which simply relate to the physical spin fluctu- 
ations of the two layers, and are coupled at the gaussian level 
via the interlayer exchange interaction. The crucial complex- 
ity of the problem comes from the fact that the two fields pos- 
sess different dynamics with different dynamical exponents, z, 
in the absence of interlayer coupling. This is a direct conse- 
quence of the fact that while one layer is open to an external 
substrate (z — 2 mode), the other remains closed from any 
external bath (z = 3 mode). As such, the bilayer system ex- 
hibits multi-scale quantum criticality. We find that, because 
the fields couple linearly, the infrared properties of the system 
are governed by a z — 3 fixed point. 

We present a perturbative renormalization group analysis of 
the system both in and out of equilibrium. In the equilibrium 
case the full system is held at a common temperature, T. Out 
of equilibrium, we introduce an effective temperature, r e jf (as 
in Ref. HI), which parametrizes the decoherence induced by 
the nonequilibrium drive. We study the flow of the system in 
the vicinity of the z — 3 equilibrium gaussian fixed point. We 
identify the line of crossover from the quantum disordered to 
the quantum critical regime in terms of the bare parameters 
of the system, and compute both temperature and nonequilib- 
rium corrections to the correlation length of the critical fluc- 
tuations in the quantum critical regime. We find that while the 
temperature corrections to the correlation length in thermal 
equilibrium gain two contributions that reflect the presence of 
the z — 2 and z — 3 dynamics, the nonequilibrium contribution 
contains only one correction that reflects z — 2 physics. This 
result signifies the fact that the nonequilibrium drive is applied 
only to the bottom layer whose corresponding fluctuations in 
the decoupled limit possess z-2 dynamics. This shows that 
the energy scale identified as effective temperature in a previ- 
ous work (Ref. does not strictly apply to this problem. 

We also consider the Langevin dynamics obeyed by the 
critical fluctuations in the quantum critical regime where the 
system effectively becomes classical and the quantum fluctua- 
tions can be integrated out. Here, the analysis is carried out in 
the eigenbasis that diagonalizes the gaussian effective action. 
Because one of the eigenmodes remains massive at the crit- 
ical point, the long-time dynamics of the full system can be 
described by a single field. We solve the Langevin equation 
for this eigenmode both in and out of equilibrium. While z-2 
physics seem to play no role in the long-time dynamics of the 
critical eigenmode in equilibrium, its effect is important in the 
nonequilibrium case. In the latter case, the Langevin dynam- 
ics display a hybrid effect of both z-2 and z — 3 physics. 
In either case, the obtained dynamics differ from any of the 
dynamical unversality classes considered in Ref. |H 

The paper is organized as follows. In Sec. [II] we begin by 
introducing the bilayer system, and present a theory to model 
the system. By integrating out all fermionic degrees of free- 
dom, the effective two-field bosonic action is derived in Sec. 
UTTl We provide a brief mean-field analysis of the effective 
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FIG. 1 : Edge-on view of the model system considered in this work. 
The two itinerant electron films (2DEGs) are coupled via short- 
ranged spin exchange interaction. The bottom 2DEG is tunnel- 
coupled to a substrate and is driven out of equilibrium by a uniform 
electric field E applied in the inplane direction. The top layer in turn 
is driven out of equilibrium by its interaction with the lower layer. 
There is no tunneling of electrons between the two 2DEGs. 



action, and establish the fixed point in Sec. [TV]for the pertur- 
bative renormalization group analysis which follows in Sec. 
IVl The correlation length is calculated in Sec. [VI] and the 
dynamics obeyed by the critical fluctuations in the quantum 
critical regime is discussed in Sec. IVHI Finally, we conclude 
in Sec. rvTTTl 



II. SYSTEM AND MODEL 

The system of interest, shown in FigQ] consists of two 2D 
itinerant electron systems which are tuned close to a spin- 
density wave instabilit y 10 n . Throughout this work, we as- 
sume the systems to be paramagnetic but are nearly unstable 
toward an itinerant ferromagnetic state with an order parame- 
ter symmetry of Ising nature. We assume no particle exchange 
between the two layers, but allow magnetic fluctuations in 
the layers to interact via short-ranged ferromagnetic spin ex- 
change. The bottom layer is tunnel-coupled to a substrate and 
is driven out of equilibrium by a uniform electric field E ap- 
plied in the inplane direction. 

The total Hamiltonian of the system is given by 



H — Hi 7 + H/ + Hb-t + H 



sub 



H, 



b—sub* 



(i) 



where b (f) labels the bottom (top) layer and sub denotes the 
substrate. The Hamiltonian for each layer has a kinetic part 
and an onsite ferromagnetic exchange interaction term with a 
common interaction strength U . They can be written as 

H b = J] £ k-^/ccl Ka .c bKo -- — Y^iSlf; (2) 



(3) 



where Cb( t )k,cr is the annihilation operator for the bottom (top) 
fermions with inplane momentum k and spin <x. The usual 
quadratic dispersion, = \k\ 2 /2m, is assumed in both layers. 
S z a i ~ Yicr a ~ c \i <T c aicr' is the magnetization for layer a, where 
i labels the inplane position, and the usual factor of 1/2 was 
absorbed in the interaction U. e is the electron charge. We use 
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a gauge in which the electromagnetic potential A = -cEf, and 
set Ti — 1 . 

Eqs.© and (0) have been written in a form which explic- 
itly breaks SU(2)-symmetry. Here, we are assuming that the 
itinerant magnets possess an easy axis so that their magnetiza- 
tion prefers to point up or down with respect to a certain crys- 
tal axis. Indeed, strong spin-orbit coupling in some itinerant 
magnets is known give rise to pronounced magnetocrystalline 



anisotrop 



17,18 



In this work, we also focus solely on the disor- 



dered side of the phase diagram where number of components 
for each of the order parameters are not expected to influence 
the results in an essential wa y 10 ' 11 . 

The Hamiltonian describing the inter-layer interaction 
reads 



(4) 



where + and - label the action for the forward and backward 
parts of the time-loop contour, respectively. The full action is 
given by 



S -S b +S t + S b _ t + S K sub + S b _ sub , 



where 

S' 



1 = J d'xclr 



.„ (id r + eA/cf 
id, + lib 



2m 



'ba- 



ll 



+ - | d 3 x{S z h ' K ) 2 , (9) 



'Oi + — + Hi 

2m 



U 



.+- | d 3 x(S z *f, 



(8) 



where J > is the strength of the ferromagnetic exchange in- 
teraction. The interaction may be mediated by magnetic fluc- 
tuations in the central paramagnetic insulator. We assume that 
these fluctuations possess long-range correlations in the out- 
of-plane direction (on the order of its thickness), but short- 
ranged correlations in the inplane direction, which justifies the 
locality of the interaction in the inplane direction. We refer the 
reader to AppendixlAlfor further details on how the paramag- 
netic fluctuations of the central insulator can give rise to an ef- 
fective spin exchange interaction between the two layers. The 
substrate Hamiltonian and the tunneling Hamiltonian, which 
describes the coupling between the bottom layer and the sub- 
strate, are given by 



H su b - ek .*: fl M-,o- ak A-.°" ; 

k,k : ,<r 

Hb-sub = £ ^ a k,L,<r C >>K<T + hc -> 
k,k : ,<r 



(5) 
(6) 



where a^.a- is the annihilation operator for the subtrate 
fermions, and ( parametrizes the tunneling strength between 
the layer and the substrate. The substrate is modeled as a 
Fermi-liquid bath with — + e(k~), where k- labels the 
out-of -plane momentum, and e(k~) is some dispersion not nec- 
essarily quadratic. The essential features of the bath are its 
density of states, p, which implies a broadening in the elec- 
tronic levels associated with the bottom layer, T = np( 2 (the 
escape rate of electrons from the bottom layer to the sub- 
strate), and its resistivity, which we assume to be very high 
relative to that of the layer so we may couple the electric field 
only to the bottom layer. 



in. EFFECTIVE KELDYSH BOSONIC ACTION 

To harness the nonequilibrium nature of the problem, 
we formulate our model using the Keldysh path integral 
formalismii. We begin by writing the Keldysh partition func- 
tion 



(7) 



SU = J j d\vs ] ;sy\ 



I' 



S* sub = I d*xdzaZ 



id, + — e(-id-) + /j. s 
2m 



(11) 



(12) 



Si-*, = -C j d 3 x (a^ =0 C + c£X, z =o) • d3) 

We are using x = (f, r), r being the inplane coordinate; z de- 
notes the out of plane direction. All fermionic fields are now 
expressed in terms of Grassmann variables. We have omitted 
position and time dependences from all fields for brevity, and 
assumed summation over the spin, k — + label a branch of the 
Keldysh contour on which the corresponding field resides. 
li, and fi s are the chemical potentials of the bottom layer, top 
layer and substrate, respectively. 



A. Free fermionic action 

The free fermionic action in Keldysh space can be obtained 
by setting U — J = 0. After performing a standard change of 
basisii, 



c T + c 

c + — c~ 



V2 

c + + c~ 



(14) 



V2 V2 
the Keldysh action for layer a is simply given by 



Vo= f ( 

Jk,cr 



C l C 2 
akcr ako~ 



(15) 



Here, k = (to, k) labels the energy-momentum three-vector 
and we have introduced the notation J k ^ = zZa- J j% ■ The 
top layer electron Green's functions correspond to those of 
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free equilibrium electrons, 

G m) = 1 
tk co - £ ± i0 + ' 

Gl = -27r/sgn(w)^-4), (16) 

where £f = -ju,. We note that the Keldysh Green's function 
is given in its zero-temperature form. 

The bottom layer is affected by both the electric field and 
the substrate, and this substantially modifies the Green's func- 
tions. Derivation of the Green's functions will be briefly out- 
lined here. More detailed derivations can be found in Ref. 
i and Ref. [H The first step is to integrate out the sub- 
strate degrees of freedom. For a frequency-independent tun- 
neling amplitude (c.f. Eq.©) the effects of the substrate 
can be parametrized by a constant level-broadening parameter 
F = np( 2 , which measures the rate at which electrons in the 
bottom layer escape into the substrate. Then the retarded and 
Keldysh self-energies which renormalize the bottom Green's 
functions are given by 

S« = -iT = (T*)*, (17) 
Ef = -2iT[l-2g(co)l (18) 

Here, g(co) is the Fermi-Dirac distribution of the substrate, 
where to is measured with respect to the chemical potential 
//j of the substrate. We note here that because of our partic- 
ular system geometry, the scattering into the bath dominates 
over electron-paramagnon scattering, which leads to consid- 
erable analytical simplifications. For instance, in Ref. the 
geometry is that of a thin strip which is coupled to a bath only 
at its boundary. In that case, electron-paramagnon scattering 
is crucial and one must calculate the electron and paramagnon 
distributions self-consistently. We also note that coupling the 
entire lower layer (and not only its boundaries) to an infinite 
bath allows one to assume that the Joule heat generated by the 
electric field is efficiently carried away into the bath. 

We now consider the effect of the nonequilibrium drive. At 
this point we specify the magnitude of the electric field that 
will be considered in this work. We will be comparing the 
electric field E to T, the escape rate of electrons from the bot- 
tom layer to the substrate, or rather to the corresponding time 
scale t = 1/2T. We introduce a quantity which will be im- 
portant in the remainder of this work: the effective tempera- 
ture scale set by the nonequilibrium field, r e ff. This quantity 
parametrizes the consequences of the nonequilibrium-induced 
decoherence. In this system, where a uniform current flows in 
the inplane direction, is given byi£ 

r eff = eE\'fT. (19) 

We will limit ourselves to weak fields: T^/Ef «: 1 and 
rT e g <k 1, where Ef is the Fermi energy of the bottom layer. 
Restricting ourselves to this regime and using the canonical 
momentum, the retarded and advanced Green's functions take 
their equilibrium forrrA 

Gff = \ . (20) 

bk co - % ± iT 



Here, & - e^-fib- The nonequilibrium Keldysh Green's func- 
tion can be obtained using the linearized quantum Boltzmann 
equatio n 15 ' 16 (QBE) for the lesser Green's function, 

eE ■ (V k + v k 3J G< = 2£ A bk - 2TG< . (21) 

With the usual parametrization, the lesser Green's function 
can be written as 

G bk = ifbkAbk, (22) 

where fbk is the distribution function, and Abk the spectral 
function. With this parametrization one obtains a QBE for 
the distribution function, 

eE • (V k + y k dj fbk = 2F[- fbk + g(co)] . (23) 

Following Ref. d, we simplify Eq.(f23]) by dropping the gradi- 
ent term in the left hand side, which is justified in the weak- 
field limit. One then is required to solve the following QBE: 

eE ■ vM bl = 2T[-f hk + g(co)] . (24) 

The Keldysh Green's function is then obtained by solving 
Eq.(f24b for the distribution function and substituting this into 
the standard formula, 

Gl k = \l-2f bk \(G R bk -G A bk ). (25) 

At zero temperature, the nonequilibrium electron distribution 
can be easily obtained and is given by 

fbk = 6(-co) + - [sgn(w) + sgn(E ■ v k )]e~^ , (26) 

where we have used the fact that g(aS) = 9(-u>) at T — (6 
being the Heaviside function). We see from Eqs. (F24l i and d26b 
that in equilibrium (E = 0), the bottom electron distribution 
reduces to the distribution of the substrate electrons signifying 
the state of thermal and chemical equilibrium between the two 
systems. We will hereafter assume for simplicity fi, = /u/, = 
Us = ii. We therefore drop the superscript t and b from the 
top and bottom dispersions and assume that both dispersions 
and frequency a> are measured with respect to the common 
chemical potential fi. 

We note here that the linear-response approach to obtaining 
the electron distribution function is justified here as long as 
the shift of the Fermi surface in momentum space (eEr) is 
much less than the Fermi wavevector. Indeed, we will find, as 
in similar previous works^, that the nonequilibrium field is 
a relevant perturbation and grows under the renormalization 
group flow. In this work, we restrict our flow up to the scale 
where eEr ~ kf so that Eqs.([25b and (|26] i remain valid. 



B. Interactions 

The interactions are given by 

«S£a = / d"x[JS^Sr + f(Sr) 2 + ftf?*) 2 ]. (27) 
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The interactions are decoupled by two one-component 
bosonic fields, rrfAx) and m*(x), which physically represent 
the magnetization in the bottom and top layers, respectively. 
The decoupling leads to 



-U 



1 



-{m K b f + J -m K b m K t 



J 



{ml + jm K t )Sl K \ + [b^t}, (28) 



where j = J/U is the interlayer exchange coupling normal- 
ized by the intralayer exchange interaction. At this point, it is 
useful to introduce a new basis for the bosonic fields: 



1 j 
j 1 



(29) 



The transformation is well-defined given j 2 + 1 . The fields 
reduce to the physical layer magnetizations when the inter- 
layer coupling / — > 0. Both bases give an equivalent descrip- 
tion of the problem but the analysis is simplified in the <f> basis 
because 4>b (<f>t) couples directly only to b (t) fermions, as can 
be seen from the third term in Eq. d28l ). Integrating out the 
fermions will not generate additional couplings between the 
bosons. In this new basis Eq. ( l28l takes the form 



^int 



-U 



1^. 



(30) 



The final Keldysh action for the interactions is obtained after 
performing a change of basis in Keldysh space. The corre- 
sponding transformation for the bosons reads 



V2 



q 



V2 



(31) 



with the analogous transformation for the conjugated fields, 
"cl" and "q" stand for classical and quantum Keldysh compo- 
nents, respectively. The resulting action then reads 



f Y—(c l c 2 



V a lk-q 



~wt rzK- 9 (^-^ (T ^)^ +c - c -]' w 

where U — U{\ — j 2 ). a and /? here are summed over b and t, 
and q = (Q, q) labels bosonic frequency and momentum. 



C. Integrating out fermions 

We now integrate out all remaining fermionic modes to ob- 
tain an effective bosonic theory in the <p basis. Here, we follow 
the standard procedure and expand the resulting action up to 
quartic order in the bosonic fields. This then yields 



(33) 



where the gaussian part reads 

^S = -J[K[(&vr 1 f»+fl* 9 ]<D 49 

-®l q (eT x m q + (b~i)}, (34) 

where e = j(Uvy l , v — m/n is the two-dimensional density of 
states at the Fermi energy and f x is a Pauli matrix in Keldysh 
space. The two-component vector = ( ^ 4> q a ) denotes 
the structure of the fields in Keldysh space. The hat above the 
polarization function, ti aq , denotes its matrix nature with the 
usual Keldysh causality structur o 14,16 , 



o K 



(35) 



Here, the retarded and Keldysh polarization functions are de- 
fined by 

K = ~ \\-G R aMq G K a k + Gl Mq Gt k \, (36) 
and = (n* )*. The quartic terms read 



n=4 



/II— t 
d\ ^ u a M(x)l 4 -"[<f> q *(x)l n . (38) 



Note that as advertised above, the quartic part of the action 
does not couple the bosonic fields. As a matter of fact, in 
the 0-basis, the only mixing appears in the quadratic part of 
the action. We now turn to the evaluation of the polarization 
functions, which will reveal the dynamics of the bosonic fields 
as well as the effects of the non"equilibrium drive on the latter. 

Integrating out soft fermionic modes is now known to yield 
non-analytic corrections to the static spin susceptibility and 
invalidate the standard Landau-Ginzburg-Wilson description 
for itinerant ferromagnets in relevant dimensions^. However, 
these non-analyticities are also known to cancel for the case 
of ordering field with Ising symmetry^. 



D. Polarization functions 

We begin with the evaluation of n f? , which is computed 
using equilibrium electron propagators of the top layer. The 
retarded component can then be obtained by analytic continu- 
ation of the standard Matsubara polarization function used for 
a clean itinerant ferromagne t 10 ' 11,21 yielding 



IK 



-i + tar- 



si 

'vflql 



(39) 



The Keldysh polarization function for the top layer can be 
straightforwardly obtained using the fluctuation-dissipation 
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theorem. At T — 0, it reads 



(40) 



The polarization functions, 11^, have been obtained 
elsewhere 8 and the detailed calculations here will be relegated 
to Appendix[B] Here we state the results. The retarded polar- 
ization function for |q| < (vft) -1 and T e $T «: 1 is given by 



n£ 9 = -i + 



qf -iT(n-v rf -q). 



(41) 



We note here that the condition |q| < {v f tY x was crucial in 
obtaining the z — 2 dynamical term above with a constant Lan- 
dau damping coefficient. In the limit where |q| > (vft)~ 1 , an 
expansion with respect to (IqlvfTT 1 becomes possible and the 
dynamical term for the bottom layer gets modified to z'Q/ Vj?|q|. 
We shall focus on excitation energies less than F and momenta 
less than T/vf where nonconservation due to escape into the 
leads is dominant. 

The parity-breaking drift term contains the drift velocity 
which reads 



eEr 

m 



(42) 



The Keldysh polarization function contains information about 
decoherence arising from the noise in the current, and is given 
by 



where 



Ll bq 



I{x) := 



-2i|ft|T 



r de 

J 2^ 



1 +/ 



1 eff 
\Q\ 



cos 6\ e *° s8 i . 



(43) 



(44) 



Absence of drift in the Keldysh polarization indicates that the 
noise does not drift in the presence of the current in the bottom 
layer. 

Once these polarization functions are inserted into Eq.(33) 
one obtains a theory of two coupled bosonic modes, one of 
which (<$>b) is characterized by a dynamical critical exponent 
Z = 2 and the other (O,) by z = 3. In the decoupled limit 
(e — > 0), the former fluctuations are non-conserved in nature 
because the bottom layer is open to the external bath. On the 
other hand, the top layer does not couple to any bath and is 
thus characterized by conserved fluctuations. 



IV. PRELIMINARY CONSIDERATIONS 

A. Mean-field analysis 

We begin with a mean-field analysis of the effective bosonic 
action (Eqs.d34"ll.(l3~8ll) in order to understand how a finite in- 
terlayer interaction and nonequilibrium drive affect the under- 
lying Stoner condition. Assuming a mean-field ferromagnetic 
state, whose solution is both static and uniform (i.e. Q. — > 



and q — * 0), we may write down the following Landau- 
Ginzburg theory for the bilayer system, 



+ Ub<p b + u 



4 

tH>t ■ 



(45) 



Here, 5 - l/(Uv) - 1; we have set the area to unity. It is 
clear from Eq.(l45b that the interlayer coupling j gives rise to a 
renormalization of the Stoner criterion. However, the (parity- 
breaking) electric field does not give an additional renormal- 
ization as the relevant corrections vanish for Q — > and 
q^O^ 

When e = 0, i.e. j = 0, the two fields decouple and the 
usual Stoner condition for the instability (i.e. l/U c = v) is 
regained for each of the two layers. However, for e > 0, the 
Stoner condition for a critical U is shifted. To see this, we 
begin with the saddle-point conditions for the theory, 



- €(/), + 2u b 4>l = 0, 

- ecpb + 2u t 4>] = 0- 



(46) 
(47) 



In general, Ub + u t , because these constants microscopically 
arise from particle-hole fluctuations, and the fermions in the 
bottom and top layers possess different dynamics. We see 
that in the mean-field equation for a given layer, the compli- 
mentary layer enters like a magnetic field that tends to align 
both order parameters. Decoupling the equations leads to 

n i xl i \_a , Si , S 2 \j3 , iS 2 »„«,-, j5 , 3<5nj«a ±1 , "lus ,a 

= (6'-/€-£)(f> a +^(u a +^Ua)4> i a +— J- JL <f>l+— ■§— C+irC- 
where a is the label complimentary to a. We see that for 5 > e 
the system is disordered and that the transition occurs at 5 = e 
as one increases U. The modified Stoner condition for the 
instability is then given by 



1 

U c = --J. 
v 



(48) 



We see that the critical U is reduced from the one correspond- 
ing to the single layer Stoner criterion because the interlayer 
interaction promotes the ferromagnetic state. The considera- 
tions above lead us to introduce the tuning parameter for the 
transition A = 5-e. The mean field phase diagram is shown in 
Fig. |2 Throughout this work, we will be considering proper- 
ties of the system in the vicinity of the transition and in the dis- 
ordered region of the phase diagram where A > 0, i.e. 5 > e. 
In addition, we will be assuming 6 > e > 0, and A «C S, e, 
which corresponds to a region close to the transition line in 
region 1 of Fig. [2] 



B. Fixed point 

The effective bosonic action (Eqs.(l34l.(l38l) consists of two 
coupled modes, <pi, and <p r , that are characterized by dif- 
ferent dynamical critical exponents. In the decoupled limit 
(e — > 0), the <p\, spin fluctuations are non-conserved in na- 
ture because the bottom layer is open to the substrate, and 
this gives rise to its z — 2 dynamics. 0, fluctuations, on the 
other hand, are conserved with z — 3 dynamics because the 
top layer is not subject to particle exchange with an exter- 
nal bath. For finite e, these modes are coupled at the gaus- 



Ordered 
(A<0) 




FIG. 2: Mean field phase diagram. The solid line separates the or- 
dered (A < 0) and disordered (A > 0) phases. The dashed line cor- 
responds to U = J; in region 1 (2), e > (< 0). In this work, we 
concentrate on region 1 of the disordered phase, near the transition 
line. 



sian level. As such, the model exhibits multi-scale quantum 
criticality. The phenomenon has been addressed in the con- 
text of the Pomeranchuk (or nematic) instability in two spatial 
dimensions 22 '-, and in the study of quantum phase transitions 
that spontaneously develop ferromagnetic order in a helical 
Fermi liquid 2 -, the latter describing the low-energy physics of 
surface states of a three dimensional topological insulator. In 
both situations, the low-energy bosonic theory consists of two 
modes, one longitudinal and the other transverse, with differ- 
ent dynamical critical exponents. The lowest order coupling 
appears at the quadratic level, i.e. <p 2 ± (p 2 , in contrast with the 
linear coupling present in our situation. In both the Pomer- 
anchuk and helical liquid scenarios, the longitudinal mode is 
Landau-damped with z = 3 and the transverse one is ballis- 
tic and undamped with z = 2. In Ref. it was argued that 
the low-energy behaviour of the system at T = is governed 
by the undamped z = 2 mode since it possesses the smaller 
effective dimension. This claim is supported by the fact that 
renormalizations due to fluctuations to the mass and the vertex 
are dominated in the infrared by the z — 2 mode. 

To carefully determine the nature of the equilibrium fixed 
point at T — 0, let us consider the renormalization to the mass 
in perturbation theory. We begin with our gaussian action in 
the Matsubara form 



FIG. 3: Tadpole diagram for the mass renormalization. The curly 
lines correspond to magnetization fluctuations <p a belonging to the 
same sector, b or t. As noted before, there's no mixing at the quartic 
level in the d> basis. 



off-diagonal mass, e, does not gain corrections from quartic 
fluctuations, the diagonal mass, 5, does gain renormalizations 
corresponding to the usual tadpole diagram of Figj3] For the 
<p b sector, this is given by 



A6 b 



= f- 

J (2x 



(27r) 2 2tt 



(50) 



The contribution of the z — 3 fluctuations near the critical 



point can be estimated by setting Q 
essentially reduces to 



|q| 3 . Then Eq 



d 5 q 



2A + |q| : 



A 3/2 , 



(51) 



where A = 5 - e is the distance to the critical point, and we 
have assumed 5 ~ e. For z — 2 fluctuations, Q. ~ |q| 2 , and the 
contribution can be estimated as 



AS b l =2 



I 



d 4 q 
2A + |q| 



(52) 



In the vicinity of the critical point (as A — > 0), the z — 3 fluc- 
tuations are expected to dominate. A similar calculation can 
be carried out for the <p, sector leading to the same conclusion. 

The key difference here from the analysis in Ref. |23| is 
that the modes are coupled at the gaussian level, and that 
the two dynamical terms enter together in all of the propa- 
gators. It is then clear that close enough to the critical point 
(where the characteristic length scale of the bosonic fluctua- 
tions £ > vj?t), the conserved dynamical term (Q/vf|q|) will 
dominate over the non-conserved counterpart (Qt), and the 
low-energy long-wavelength theory will be governed by z — 3 
dynamics. We now study the renormalization group flow in 
the vicinity of this z — 3 fixed point. 



c(2) _ 



fbq <t>t q ) 



-^bq 



4>ln 



(49) 



where £^ q 



5 + \q\ z + |0|T and £^ = 5 + |q| 2 + |0|/v F |q| 
are the equilibrium inverse susceptibilities. The propagators 
for the modes, obtained by inverting the 2x2 inverse sus- 
ceptibility matrix, are (<pb(q)<f>* b (q)) = -C^/D^, i4>MWM)) = 



and {4> b {q)f t (q)) = {<f>,(q)f b (q)) = e/D^, where 
= J^X^ - e 2 is the determinant of the matrix. While the 



V. RENORMALIZATION GROUP ANALYSIS 

To aid with the analysis we begin with the effective gaussian 
action (c.f. Eq.d34b). 



j( 2) 
'eff 



-er x L t j \ ®tq 
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where X. a = (Uv) l t x + fl aq . The explicit expressions for the 
matrix elements read 



Ll = 6 + |q| 2 + i\ d ■ qr - iQ.T, 



-2i|Q|r 



and 



Xf = 5+|q| 2 -/ 



Q 
Vflql' 



Lf = -2i 



|Q| 
Vflql 



(54) 
(55) 

(56) 
(57) 



The quartic action was given in Eq.((38). 



A. Flow equations out of equilibrium 



We study the flow in the vicinity of the z — 3 gaussian 
fixed point with T e s = A = u" — 0. Here, T — through- 
out. At this fixed point the most relevant dynamical term (i.e. 
Q/Vflql) remains marginal. In Eq.d53b. we have defined our 
momentum cutoff A, which in turn defines an associated cut- 
off energy scale v F A. All energy quantities will be considered 
to be in units of this scale. The bare parameters are taken 
to be A, T sS , u", 8, e <K 1 with 5 ~ e » A. We begin by 
integrating out fluctuations whose modes reside in the shell 
A./b < q < A. Here, b > 1 is the scaling variable. In order 
to preserve Keldysh causality, we integrate over all frequen- 
cies (i.e. -oo < O < oo) at every mode elimination step. 
We next rescale momentum, q — > q' lb, to restore the cutoff 
back to A. At the z — 3 fixed point, the frequency scales as 
Q — » Q'/fr 3 and this keeps the coefficients of |q| 2 and Q/vj?|q| 
terms invariant under the flow provided that the fields scale 
as — > b 1 ^[(ffaqY ■ The (r e ff/£2)-scaling present in the 
Keldysh components suggests that r e jj is relevant at the fixed 
point and scales with dimension of 3, i.e. — > TL./b 3 . The 
non-conserved dynamical term is irrelevant at the fixed point, 
and this is reflected in the scaling r — > r'b. Then, the re- 
quirement to maintain the drift term invariant under the flow 
dictates the scaling of the drift velocity as Vd — » v'Jb 2 . Note 
that the drift term contains the combination tvj, which has 
scaling dimension 1, just as in Ref. HI Both the diagonal (6) 
and the off-diagonal (e) mass terms scale with the usual di- 
mension of 2. Finally, at the z = 3 fixed point, all quartic 
coupling constants are irrelevant with [u"\ = — 1. Therefore, 
in the standard lor o the order parameter self-interactions 
are not expected to play a significant role in the effective low- 
energy theory. 



RG equations to lowest (linear) order in the quartic coupling, 

dT eS (b) 



dlnb 

de(b) 

dlnb 

dv d (b) 

dlnb 

dr(b) 

dlnb 

du a „{b) 

dlnb 

d5 a (b) 

dlnb 



3T eff {b), (58) 

2e{b), (59) 

2v d (b), (60) 

-T(b), (61) 

-<0>) + 0([<] 2 ), (62) 

26 a (b) + 3u a l (b)f a (.T eS (b), r(b)). (63) 



We will distinguish renormalized parameters from the bare 
ones by explicitly writing the ^-dependence for the former. 
In Eq.d63l, a labels b and t. Indeed, the bare values satisfy 
6i, = S, — 6; the label denotes the fact that these mass terms are 
subject to different renormalizations from fluctuations. The 
function, f a (T e g(b), r(b)) is given by 



f a (T eS (b),T(b)) = iK : 



I 



da 

27 



^d,o), 



(64) 



where K2 — J j^rdiq- 1) = 1 /2n, we have set the momentum 



cutoff A = 1, and 



r K If* I 2 + £ 2 f K 

ap ~ I rR rR _ e 2\2 
^bp-^tp fc 1 



(65) 



is the Keldysh Green's function for the mode <p a , it was 
obtained by inverting the inverse susceptibility matrix of 
Eq. (l53l l. Solving these equations, we arrive at 



T eB (b) 
e(b) 
v d (b) 

T(b) 

u° n (b) 



= Teffb* 

= eb\ 

= v d b 2 , 

= Tb- [ , 

= uZb~\ 



(66) 
(67) 
(68) 
(69) 
(70) 



8 a (b) = b 2 [6 + 3 



r*\n b 

Jo 



dxe- 3x f a (T eS e ix ,Te- x )].ai) 



B. Flow equations in equilibrium 

To elucidate the nonequilibrium problem at later stages, we 
also consider flow equations for an equilibrium situation in 
which the entire system is at temperature T. Here, the electric 
field is set to zero, thus r e ff = 0. In this case, the gaussian 
action, Eq.(|53"ll, becomes 



Z R b = 6 + \q\ 2 -i£lT, 



To summarize our scaling analysis, we write down the set of 



Ll = -2/Qrcoth 



(it)' 



(72) 
(73) 
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and 

J* = 5 + \q\ 2 -i-^-, (74) 

v F \q\ 

£f = ~2i— -coth — . (75) 
Vf|q| \2T) 

An analogous renormalization group analysis as in the 
nonequilibrium case can be applied here, and the solutions 
to the corresponding flow equations become 



Tib) = Tb s (76) 

e(b) = eb\ (77) 

r(b) = rb-\ (78) 

u a „(b) = u a n b~\ (79) 



rAn b 

6 a {b) = b 2 [6 + 3u1 dxe- 3x f a (Te 3x ,Te- x )]. (80) 
Jo 

The /"„ -functions are still given formally by Eqs.(f64b.(l65l> with 
the replacement r e ff — > T, and the equilibrium inverse suscep- 
tibilities, Eqs.(|72]i-(|75]l, must be used. 



VI. CORRELATION LENGTH 



A. In equilibrium 

We begin by considering the correlation length in equilib- 
rium. We define the quantum-to-classical crossover line for 
the system, and compute the leading temperature correction 
to the correlation length in the quantum critical regime. 



1. Quantum disordered regime 

In this regime, the /^-functions are computed at T = 0. 
This gives a temperature independent shift proportional to u" 
to the two diagonal masses. We then have 

Jr»oo 
dxe- 3x f a (0,Te- x )], 
o 

:= b 2 (6 + A6° a ). (85) 

The correlation length in this regime is then given by 

ST 1 * A + AS° b + A<5? =: r, (86) 

where we have computed to lowest order in A/ 6 <k 1 and 
A6®/5 <sc 1. To find the condition on T for the occurrence of 
the regime, we impose that A(b') = 1 while T{b') < 1. Recall 
that A(b) « A bh (b) « A„(b) ~ A h ,(b). This then translates to a 
condition on the bare quantities, 

T < r V2 . (87) 



In the current bosonic basis, there are three two point cor- 
relation functions one can study, each with its own associated 
mass scale, which we define as 

Kp = <^ ? >l ? =o> (81) 

where q is a frequency-momentum 3-vector. In the above 
equation we have omitted the Keldysh indices because one 
can use an equilibrium formalism to define the mass scales as 
they don't depend on the nonequilibrium drive. Note also that 
At, t = A,/,. Within the gaussian theory, the bare mass scales 
read 

5 2 -e 2 

Abb = A„ = — — , (82) 
o 

6 2 -e 2 

A b , = . (83) 

e 

In the relevant limit, 5 > 0, e > 0, and 5 > e » A, 
Eqs.(|82]i,([83]l become 

A bb = A„ * A bt * 2A. (84) 

Here, the mass scales were computed to lowest order in A/ 6 <K 
1 and A/e <k 1. One sees that they all vanish as the transition 
is approached (A — > 0). This allows one to then define the 
physical correlation length of the system, via £,~ 2 := A. We 
now compute £ both in and out of equilibrium focusing on the 
corrections arising from either T or T e ff . 



2. Quantum critical regime 

When the inequality in Eq.(|87b is reversed, the characteris- 
tic energy scale of the fluctuations, 1 : , become smaller than 
temperature and the fluctuations become classical. Here, we 
compute the leading temperature correction to the correlation 
length at scale b", where T{b*) = 1. Since 5 ~ e » r ~ T 2/3 , 
both 6 and e become much larger than 1 at some stage of the 
flow. We define the scale b\ as 5{b\) ~ e{b\) = 1. Then the 
integral in Eq.(f80b must be split into two regions: 1 < b <b\, 
and the other b\ < b < b*. In the former region, the integral 
can be computed assuming 5 ~ e <K 1 inside Eq.(|65l). In the 
latter region, the integral is computed assuming 8 ~ e » 1. 
The detailed derivation will be relegated to AppendixICl here 
we state the results. We define 

A5 a := 3< f dxe~ 3x [f a (Te 3x ,Te~ x ) — f a (0,re~ x )]. 
Jo 

(88) 

Here, one integrates up to b* by splitting the integral as adver- 
tised. One then arrives at the following temperature correction 
to leading order, 

M ° = £H 1 + !^ /3 )- (89) 
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The correlations length in the quantum critical regime is then 
given by 

u b + u\ I 3t m \ 

^ Kr + ^2 ±T [ l + Y T )■ (90) 

The linear-temperature correction in Eq.d90b is consistent with 
the corresponding correction obtained in Ref. [Tl| for the case 
d — 2 and z = 3. The T 4 ^ 3 -correction, which is proportional to 
r, is an additional correction that arises because of the z — 2 
dynamics present at the z — 3 fixed point. Both the leading 
z — 3 correction and a sub-leading z = 2 correction enter into 
the correction because both modes are subject to a common 
temperature T. 



B. Out of equilibrium 



perature T can as shown in Eq.(|9Q|>. However, the lack of a 
linear- T e ff correction for the correlation length indicates that 
r e jf cannot be identified strictly as an effective temperature of 
the system. 

We note that in equilibrium two corrections that reflect the 
presence of z — 2 and z = 3 dynamics enter into the corre- 
lation length. The dominant contribution (i.e. linear- T con- 
tribution) arises essentially from the z — 3 dynamical terms 
in the regime of interest. In systems with multi-scale quan- 
tum criticality, modes with different dynamics often nontriv- 
ially contribute to various thermodynamic quantities in differ- 
ent bare parameter regimes. In thermal equilibrium, interest- 
ing crossover behaviour in thermodynamic quantities in a sys- 
tem with multi-scale quantum criticality has been addressed 
in the context of Pomeranchuk instability 2 - and the Kondo- 
Heisenberg model 2 ^. 



We now compute the leading nonequilibrium correction to 
the correlation length at scale b*, where T e ${b*) = 1. We now 
assume T — 0. Here, we only concentrate on the quantum 
critical regime because the correlation length in the quantum 
disordered regime is still given by Eq.(f86b. Of course, the 
condition for the occurrence of this regime now reads T e s < 
r 3/2 . 

The leading correction due to the electric field can be com- 
puted by splitting the scaling regimes into 1 < b < b\ and 
b\ < b < V as in the equilibrium case. We then obtain, 



3u l l 
An 3 



A6 a « ^\rT%\ 



(91) 



The correlations length in the quantum critical regime is then 
given by 



neq 



r + 3- 



u b + u\ 



4tt 3 



< 3 - 



(92) 



Once again, the detailed calculations are presented in Ap- 
pendix [C] In contrast to the equilibrium result, Eq.(l90t. the 
correlation length gains a correction of order tB only. The 
difference arises because in the decoupled limit the two lay- 
ers are at different effective temperatures. This can be read- 
ily checked by comparing the two Keldysh polarization func- 
tions, Eqs.d40l) and d43l . One sees that the top layer is at zero 
effective temperature while the bottom possesses a finite T e ff 
because the field directly couples to the latter. Once the in- 
terlayer coupling is restored both <pt, and <p, sectors feel the 
effect of the field as it is evident from Eq.(|92l. However, due 
to the absence of an effective temperature in the top Keldysh 
polarization function, a linear- correction does not arise in 

Results in this section suggest certain similarities between 
the role of T e $ in the nonequilibrium case and temperature T 
in the equilibrium case. Both T e $ and T induce a quantum-to- 
classical crossover, across which the behaviour of critical fluc- 
tuations evolves from quantum to effectively classical. The 
crossover energy scale is defined in Eq.(|87b where T can also 
be replaced by 7^. Furthermore, Eq.(|92l implies that r e fj 
can cutoff the divergence in the correlation length just as tem- 



C. Eigenmode analysis 

We have found that one must use caution in determining the 
mass scale which relates to the inverse square of the physical 
correlation length, for neither 5 nor e alone defines this scale. 
A sensible procedure would be to determine the eigenmodes 
of the system that diagonalize the gaussian action Eq.(t34b. and 
directly read off the corresponding masses by observing the 
eigenvalues of the action. This analysis was carried out, and 
is presented in Appendix ITJ1 We find that while one of the 
eigenmodes becomes critical at 5 — e, the other remains mas- 
sive. We then integrate out the massive mode and construct an 
effective critical theory in terms of a single critical mode, A q . 
In the appendix, we carry out a renormalization group pro- 
cedure on this single-component critical theory and the cor- 
responding T and corrections to the correlation length of 
this field is computed. We find that the correlation length ob- 
tained in the eigenbasis gives identical results to those pre- 
sented in Eqs.d90l.d92l 



VII. ORDER PARAMETER DYNAMICS 

We now consider the dynamics obeyed by the critical fluc- 
tuations in the quantum critical regime where the theory is 
effectively classical and the quantum fluctuations can be in- 
tegrated out from the effective low-energy theory. The stan- 
dard procedure to obtain the dynamics is outlined in Ref. 
and has been applied in RefsOH in establishing nonequilib- 
rium dynamical universality classes. In a completely different 
setting but in a related work, order parameter dynamics near 
phase transitions in a driven interacting bilayer lattice gas has 
also been addressed in RefsElzzl 

We begin by decoupling the terms quadratic in the quan- 
tum Keldysh fields by making use of a second Hubbard- 
Stratonovich transformation. The decoupling field plays the 
role of a noise source for the classical fluctuations. Finally, we 
integrate out the quantum fluctuations and obtain a stochastic 
equation, or Langevin equation, for the classical fluctuations. 
We hereafter drop the "cl" subscript from the fields. All cou- 
plings will also be assumed to take their renormalized values. 
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In the <p a basis, the system of coupled, linear Langevin 
equations is 



-CL ~ e \( 4>bq \ _( TJbq 
\<Ptq) \ TJtq 



-e L R 



(93) 



where r\ aq are the Fourier components of the noise sources; 
their correlators are given by the Keldysh inverse susceptibil- 
ities 

(VaqTlPq' > = <V<% + l'Waq ■ ( 94 ) 

One can obtain the noise-averaged correlators for the fluctua- 
tions by inverting the coefficient matrix: 

\£ R I 2 -2 

^tql K 



<^>=l^nf 9 + — (95) 



where D q = -£*.-£^ - e 2 is the determinant of the coefficient 
matrix. We mention here that for finite interlayer coupling 
all bosonic correlators feel the presence of the parity-breaking 
drift term because they all contain J3. One can repeat the 
same procedure for the critical eigenmode (see AppendixlDb. 
whose Langevin equation reads 



(98) 



with the noise correlator y]qTlq'/ = V(k^(9 + <?')■ From Ap- 
pendix|D] we have that 



Xk 1 = \(Kl q + U%). 



(99) 
(100) 



The noise-averaged correlator for the eigenfield then reads 



nf ? ) 



kU R bq +Zfq)-e\- 



(101) 



We now argue that the three correlators in the basis and the 
one for the critical eigenmode contain the same physics in the 
regime of interest. Indeed, we are considering A as 5 - e «c 
5, e, such that at the end of scaling 5, e > I. Thus we can write 
\£aq\ Z ~ Neglecting differences between 5b and 5, due 
to one-loop renormalizations, we have ((paqfip) ~ (AqA.*)/2. 
Hence, we will study the fluctuation dynamics in terms of the 
critical eigenmode, A. 

In equilibrium, Eq.d98l) in momentum-time space becomes 



1 



Vflql 



3,A(q,f) = -(A + |q| 2 M(q,f) + 77(q,f), 



(102) 



where we have dropped the z — 2 term from the time- 
derivative because we are in the regime of |q| < F/v/r. In 



the quantum critical regime, the noise correlators become 

IT 

<77(q, fMq', f')> * —r-M - f)6(q + q'). (103) 
vplql 

We then find the noise-averaged correlation function for the 
fluctuations 



{A(q,t)A(-q,t')) * - e -^'-'\ (104) 

In the out of equilibrium case, the Langevin equation be- 
comes 



1 



v F |q 



d,A(q, t) = -(a + |q| 2 + l -\ d ■ qrj A(q, t) + T](q, t). 

(105) 



In the limit of large T e ff , we find 
2T ef fT 



<77(q,f)77(q',f')> 



6(t-t')6(q + q'). (106) 



The noise-averaged correlation function then reads 



(A(q,t)A(-q,t')) 



T e gTv F \q\ _, 

7lA S 



f (q|(A+^v d -qT)|f-?'| 



(107) 



In equilibrium, the correlation function describes fluctua- 
tions which are conserved in nature with the decay rate van- 
ishing in the long-wavelength limit as ~ |q|. Because the ef- 
fects of temperature are felt equally by both bosonic modes, 
the z = 2 dynamics play no role in the long-wavelength the- 
ory of the order parameter dynamics. The dynamics obeyed 
here is essentially identical to Model B in Ref. |^ but with 
the Landau damping parameter replaced by Vf|V|. The situa- 
tion is different in the nonequilibrium case. Here, the delta- 
correlated noise contribution only arises from that reflects 
z = 2 dynamics since the driving field (the source of noise) is 
only coupled to the bottom electrons. The Langevin equation 
then displays a hybrid effect of both z — 2 and z = 3 physics; 
while the white noise correlator reflects z — 2 physics, the 
damping is governed by the z — 3 dynamical term. As a result, 
the corresponding noise-averaged correlation function for the 
fluctuations yields a unique behaviour where long-wavelength 
correlations vanish as ~ |q|. In addition, the advertised drift of 
the flucutations can be explicitly seen in Eq. J 107b . 



VIII. CONCLUSION 

In summary, we considered critical properties of a nonequi- 
librium bilayer system of itinerant electron magnets. Starting 
from a microscopic fermionic model subject to an external 
drive, we derived a coupled theory in terms of two bosonic 
fields which are related to the physical magnetization fluctua- 
tions of the two layers. In the limit of no interlayer coupling, 
the fields obey different dynamics with different dynamical 
critical exponents (z = 2 and z = 3), leading to multi-scale 
quantum criticality in the coupled system. We found that the 
applied current leads to both a drift of the magnetic fluctua- 
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tions in the coupled bilayer system and to decoherence. The 
latter phenomenon is more subtle in our system compared 
with the analogous single layer scenario because the nonequi- 
librium drive is applied to only one of the layers. This causes 
it to play a role distinct from temperature in the thermal equi- 
librium case where both layers are held at a common tem- 
perature. The differences are illustrated by comparing tem- 
perature and nonequilibrium effects on the correlation length 
and the dynamics obeyed by the critical fluctuations. A cru- 
cial feature of the work is that the two fields couple linearly. 
We found that the infrared properties of the system then are 
governed by the dynamics corresponding to the higher effec- 
tive dimension. In light of Ref. |23[ it would interesting to 
consider a coupled order parameter theory where the fields 
couple quadratically in the context of itinerant electron mag- 
nets. In this case, the effective theory possesses a discrete 
Z2 x Z2 symmetry corresponding to independent transforma- 



tions m a — > — m a for a — b,t ensuring that no linear coupling 
can be generated during the RG transformation. If both fields 
become critical simultaneously, the low-energy properties are 
then expected to be governed by the field with the lower ef- 
fective dimension. 
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Appendix A: Interlayer spin exchange 

We now consider how the central insulator can mediate an effective spin exchange interaction between the top and bottom 
layers. Imagine we have two 2D itinerant ferromagnets sandwiching a thin but 3D insulator with thickness L. For concise- 
ness, we envisage the insulator as a quantum Ising paramagnet in the vicinity of a quantum phase transition to a long-ranged 
magnetically-ordered phase. The continuum quantum field theory for this model can be expressed in terms of a coarse-grained 
order parameter fluctuation field, <p, subject to a confining potential. For simplicity, we truncate this potential at quadratic order 
with some invertible dynamic correlation matrix for the fluctuations, xfc< Xj, t), where x, = (r,-, z,) labels the lattice sites in the 
3D insulator. The QFT in the vicinity of the magnetic quantum critical point is then given by 28 

Sins = I dT^<t>(Xi,T)x~ l (Xi,Xj,T)<l>(Xj,T). (Al) 
J U 

Since we are only interested in how this central insulator generates an effective spin exchange, we will work in the equilibrium 
formalism for simplicity. 

We now consider the interaction terms between the insulator and the two itinerant magnetic layers. We assume that the plane 
of the bottom layer is z = and the top layer is at z — L. The simplest plausible interaction would be a ferromagnetic onsite 
exchange (K > 0) between the layers and the insulator: 

St-bm = -K j dTj^S*(ri,T)(f>(Xi,T)d ZirL , (A2) 
S h - ms ^-KjdT^ §l(r h T)4>(Xi, T)6 Zh0 . (A3) 

Here, S] and S z b are the Ising spin fluctuations in the top and bottom layers, respectively. We integrate out the central fluctuations 
tp assuming that their in-plane correlations are short-ranged to obtain the following effective action 

Seff= f rfr2[-^;rt£.W(ri,T)S^ (A4) 

r, 

where we have assumed that the inplane correlations of the insulator are short ranged and independent of r. x( x i> x j> T ) = 
S Th TjX{Zi,Zj). As advertised, an effective interlayer spin coupling J = K 2 [x(L, 0) + x(0,L)] is generated. We see that -K 2 x(L,L) 
and -K 2 x(0, 0) lead to finite renormalizations of the intralayer exchange term. 
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Appendix B: Polarization functions for the <p b sector 

In equilibrium, the retarded polarization function in the (pi, sector is given by 

n£ * -1 + Iql 2 - iOr, (Bl) 

where we are assuming Qr < 1 andqTVf < 1. In the regime where qrx'f » 1, the dynamical term becomes Q./vfq, i.e. conserved 
in nature. This is because on sufficiently short length scales the electrons are unaffected by the presence of the substrate. The 
equilibrium contribution to the Keldysh polarization is given simply through the fluctuation-dissipation theorem, 

n£ = -2»m|T. (B2) 

To get the nonequilibrium corrections to the polarization functions we use the nonequilibrium Green's functions for the bottom 
electrons obtained in Sec lHI Al From Eqs.d20b.d25l>. the nonequilibrium contribution to the retarded polarization reads 



an«(P,io*-i— f -r^ 



eE ■ VkT 



(4k + r2 ^k + r2 )' 
Expanding this result for small q one finds 



The nonequilibrium contribution to the Keldysh polarization is given by 

r -i C -2iT -2iT ( a_ 

OTf(O,0) = - - ,, 2 ^r2 7 TTj-pa -sgn(w + n)sgn(w)« ^ 



(B3) 



OTZ(P,q)"i-Li«E.q. (B4) 



- sgn(oi) sgn(ai + fi)e l' E - v H + [sgn(w) sgn(w + O) + l]e W^R e | . (B5) 

Here, we have linearized the fermion spectrum in the exponents. For a weak field (i.e. <K r _1 ) the functions in {. . . } in 
Eq. (IB5l ) are strongly peaked at a> values where the numerators in the exponents vanish. Making use of this we obtain 

4/p2 r i 

ffllf (fl, 0) * — 5 — {. . . }, (B6) 

where {. . . } is unchanged from Eq.( |B5t . This makes the w-integral trivial, yielding the result 

" 2 C kdkd9T etf \cos9\e T d*»* 

where 9 is the angle subtended by the Fermi velocity and the electric field. Performing the momentum integral gives 



k „ 4r C kdkd9T ef f\co&9\e T **»<* 
OTf £2,0 =-/— - ett ' 2 ' , (B7) 



«lf(O,0)«-2jT|O|/|^J. (B8) 



Appendix C: Evaluation of the correlation length 

In this appendix, we present the derivation for the leading correction to the correlation length in the quantum critical regime. 
Here, the analysis is conducted in the (pb~<Pt basis. The results are summarized in Sees. IVI A2l and lVIBI 

Recall that 5 a (b) in equilibrium is given by Eq.(l80b. The correction to the zero-temperature distance from criticality, r, is then 
given by 

k8 a = 3ul[ dxe- 3x [f a (Te 3x ,Te- x )-f a (Q,Te- x )l (CI) 
Jo 
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If we define Af a := f a (Te , re x ) - f a (0, re we have 



Af a = iK 2 J — ADfCT.T), 



where K 2 = l/2n and 



ADf(l,£J) 



-2/fi 



ADf(l,Q) = -2ifi 



T[(<5+l) 2 + fi 2 ] + e 2 
I (5 + 1 - iQ.T){5 + 1 - in) - e 2 | 

(£+ l) 2 + (nr) 2 + Tf 2 
\(6 + 1 - ifiT)(<S + 1 - in) - e 2 | 



(C2) 

(C3) 
(C4) 



The fi-integral gets most of the contribution from |fi| < 2T. So we can cutoff the integral, and approximate coth(fi/2!T) * 2T/CI. 
Then we obtain 



ADf(l.n) * -2i(27-|n|) 

ADf(i.n) * -2/(2r-|n|) 



r[(5+l) 2 + fi 2 ] + e 2 
\(S + 1 - iQ.T)(6 + 1 - in) - e 2 | 2 ' 

(6 + l) 2 + (nr) 2 + re 2 
\(S + 1 - /fir)(c5 + 1 - in) - e 2 P 



(C5) 
(C6) 



To scale up to V ', where r(Z?*) = 1, we split the integral into two regimes: 1 < b < b\, and b\ <b < b*, where b\ — 6 l ^ 2 . Up to 
scale b\, S, e <K 1 and T <K 1. So in this regime, we can use the following expressions 



AD* 



ADf 



-2i(2r-|n|) 

-2i(2T-|n|) 



1 + (nr) 2 ' 

1 



1 +n 2 ' 

With this approximation, the fi-integral can be done trivially. Expanding the result for small T, we get 

A/ 6 = — tT 2 +0{T\ 

71 

Af, = — -T 2 +0(T 3 ). 

7T 



(C7) 
(C8) 

(C9) 
(CIO) 



To evaluate the integral up to b\, we first do a change of variables, x — > 7\Z?) =: y, which makes the measure dx = dy/3y, and 
the lower and the upper limits T and T/S 3 ^ 2 , respectively. Then the corrections are given by 



Jt 3y y 



l/3\ 



We then obtain 



A<? 



(i) 



o 7 2/3 

7r~ o 



A^ u = — u[T 



T 2/3 



(Cll) 

(C12) 
(C13) 



To scale between b\ and b* , we use Eqs. dC5b . dC6t in the limit of S ~ e » 1. Since A «c 1 throughout the flow, we note that 
6 2 - e 2 ~ SA < 1 during the flow, and we drop this term from the denominator. Then Eqs.( |C51 l, dC6l > can be approximated by 

ADf(l,fi) * -2/(2r-|fi|) Jf 2 ^" 2]+e2 ; (C14) 

4£ 2 + (6£l(l + t)) 2 

K 5 2 + (fir) 2 + re 2 

ADf(l,n) « _2i(^_^_^__. (CIS) 
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We then find that 



t5 2 + e 2 - 1 + To 



Afb 
A/, 



4n 2 6 2 
re 2 +5 2 



An 2 



-i 1 + T t 
y^Z y^Z 



4tt 2 <5 2 An 2 

where we have used the fact that 8 ~ e. The corrections from this region of the integral is then given by 



J T/S 3,2 3y y 



T 

r, 7 



im 



We then obtain 



n° T 

A6 (2) = ^~ 
a 4n 2 



i + 3 -1t 



1/3 



Since T 2 ^ <k 6 in this regime, Eqs. (|C12t , dC13l l are sub-leading to Eq. dC19l ). We therefore obtain to leading order 

3t 



AS a 



4k 2 



l + -T ] > 3 
2 



(C16) 
(C17) 

(C18) 

(C19) 

(C20) 



In the nonequilibrium case, recall that 5 a (b) is given by Eq.dTTli. The correction to the zero-temperature distance from 
criticality, r, is then given by 



-In b' 



A5 a = 3< f dxe^ x Ua (T eS e 3x , re'*) - f a (0, re"*)] . 
Jo 

In this nonequilibrium case, we obtain 



\ l £2 l / \(6 + 1 - iQ.T)(6 + 1 - z'Q) - e 2 \~ 

9 

re' 



ADf(l,Q) = -2/|Q|7 



/ |(£ + 1 - ;Ot)((5 + 1 - in) - e 2 \ 
In the above, the drift term was dropped since it does not reach cutoff scale at b*& 



(C21) 

(C22) 
(C23) 



We first scale up to b\. As in the equilibrium case, we have 5 <K 1 and e <k 1 in this region, and Eqs.( lC22l i.( lC23t can be 
approximated by 



AD* 
AD? 



-4/rr eff l_ 

7i 1 + (Or) 2 ' 
-4/Tr eff e 2 

71 [1 +(Qr) 2 ](l +Q 2 )' 



Once again, we have used the fact that AD* is appreciable only for \Q\ < 2T e ff and that in this region, 



I mi/ mi n 



Performing the f2-integral up to 2T e ff, we then obtain to lowest order in T e g, 

Af, = ^Tl^OiTl,). 

71 



(C24) 
(C25) 

(C26) 

(C27) 
(C28) 
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The corrections up to b\ can then be computed, giving 

3 



, r 2/3 
a J £■ t4/3 eff 



3 r 4/3 

a J r 2-T-2/3 eff J f T 2 /pin\ 

' = 2^" lTe ^ ~ * 2^2«i r:r eff- (C30) 



To scale in the region b\ <b < b*, we use 



ADfd.fl) « ~ 4 ' Yreff , T<? (C31) 

h ' n 4-S 2 + [<K2(1 +r)] 2 

ADfd.fl) * ~ AlTTeS A ^ ^' „ 2 . (C32) 
t ¥+ [£Q(1 + r)] 2 



Performing the Q-integral and expanding for small T e g, we obtain 



Af b * A/, * (C33) 



where we have again used the fact that 5 ~ e. The corrections are then given by 



3™° 

A<e - -jJ-Tg. (C34) 



As in the equilibrium case, the contributions from the first region of the integral (Eqs. dC29] l,( IC30l )) are sub-leading to those from 
the second region (Eq. dC34b ). Therefore, to leading order, we obtain the following correction: 

3TU a 

A^^eT- (C35) 



Appendix D: Eigenmode analysis 



1. Determining the eigenmodes 



In this appendix, we develop a critical theory in the basis of the eigenmodes. The renormalization group analysis will be 
carried out in this basis, and the correlation length calculation will be reconsidered. We determine the eigenmodes of the system 
by diagonalizing the gaussian action Eq .(13411. What we find in the following is that the diagonalization procedure gives rise to 
two new bosonic modes which possess different masses. While one becomes critical at 6 = e, the other mode remains gapped. 
This allows one to integrate out the latter and arrive at an effective low-energy theory in terms of one critical mode. This is in 
stark constrast with Ref. |23] where both bosonic modes simultaneously become critical. 

To aid with the eventual block-diagonalization of the Keldysh action, we first diagonalize the equilibrium action in the Mat- 
subara formalism, then map the theory back to Keldysh space. The Matsubara gaussian action reads, 

where = 6 + |q| 2 + |Q|t and £, e rq = S + |q| 2 + |Q|/i'f |q| are the equilibrium inverse susceptibilities. Expanding to lowest 
order with respect to the fluctuating parts of the polarization functions, the eigenvalues of the action read ji-C^ + ± e. The 
expansion is well-defined in the interlayer coupling regime considered here (i.e. d > e > 0). If we now denote the eigenmodes 
by A\ q and A^, the gaussian action then reduces to 

(2) _ r * q da, J Ul] + £*>ii -e o \u lq \ 

^"JW^" 2 " ) { Ul + £*)I2 + e)\l 2q y (D2) 

We can now map this Matsubara action back into the Keldysh form by simply replacing both fields by two-component Keldysh 
fields, i.e. A, q — > A, q = ( Af q A^ ^ , and replacing each polarization function by its corresponding matrix. The Keldysh gaussian 
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action is then given by 



iS 



ID 
eft 



J {In) 2 2n \ u > 2 « ' 



(Z bg + X.,n)/2 - €T X 











+ £ t? )/2 + ef x j \ A 2? 



A 



i 'i 



(D3) 



where £ ff9 = (Uv) l r x + I\ aq - We now consider the quartic terms. Recall that in the tp basis, the quartic interactions were given 
by 



iS » = J <i 3 * {-/ [«? (^) 3 4>l + ^ (^) 3 ] + [«S (tf )' (tf + u\ 



+ (ft <h> t), 



(D4) 



where we have written terms in the form in which u" are all real. To transform the quartic terms to the eigenbasis, we use 
the transformation matrix in the limit of q — > since momentum and frequency dependent parts generate terms that are (RG) 
irrelevant compared with the leading constant coefficients. We then obtain 



J_/l 1 
V2 I 1 "I 



cl,q 
cl.q 



(D5) 



Inserting this transformation into Eq. (ID4b gives the quartic terms in the eigenbasis. 



The mass terms of the two eigenmodes in Eq. (ID2b are 6 - e and 5 + e for Ai and A2, respectively. This implies that while 
Ai becomes critical at 6 = e, A2 remains gapped. We may therefore integrate out the gapped mode from the theory and obtain 
a single-component effective action only in terms of the critical eigenfield. In the process of integrating out A2, we only retain 
terms that are quadratic in the latter. Terms linear and quadratic in A2 emerge from the quartic terms, however the former will 
generate contributions beyond quartic order in the critical field. After these simplifications, the quartic action becomes 



iS 



(4) 
elf 



= -i j d 3 x ui U1) 3 A* + Mi (4) 3 + j d 3 x u 2 (A? f (4) 2 + K 4 (4)' 



2?' 



(D6) 



with Qn = 3u\Af k A^ _ k _ , - iu2A^ k+q A\ _ k _ q ,- Here, m,- = (w* + m')/4. Q12, Q21 and Q22 are also terms quadratic in A\, 
however, we do not write them explicitly here because they will not be necessary in the following discussion. Now performing 
the gaussian integral over A2, and expanding the resulting Trln to linear-order in S,, we obtain 



iS fs = ~i f^l + A)/2 - rfj Ai, - i f\\ q ( ° * ) A l9> 



(D7) 



where 



and 



3j Y \ k iD 2k> m^aif iD*, 



2q 2q 
2q 







[Ubq + tq)/2 + €T X \ 



(D8) 



(D9) 



Both 771 and 772 are real quantities, and we have used the fact that D^t, t) = D^ q {t, t) = 0. We find that 112 — > as T — » so 772 
does not change the form of the Keldysh term for Ai. We will omit these terms. The u\ term gives a small renormalization to 
the mass of the critical mode. We absorb the 771 into the new mass: A = 6 - e + 771 . The full effective action for the critical mode 
is now given by 



l<Jeft 



IMxr ? $Y lq ' i \ ^(^f 4 + ^fM 3 ] + / d'x[u 2 (Atf(Afj 2 + u 4 (A^ 



(D10) 
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In equilibrium, 



Xr x =A + |q| 2 ~(nT+-^-|, 



2 \ v F \q 
Q 
v F \q 



X K l = -/coth^)^r+-^| 



and out of equilibrium, 



Xl 1 =A + |q| z + -v £/ -qr--lQT+^-|, 



Xk = ~i 



2 ' 2\ v F |q 

mi 



Vflql 



|0|r(l + /(^)). (D14) 



Note that /(T eff /P) ~ r eff /|Q| for r eff > and coth(Q/2r) ~ T/Q, for T > |Q|. We then see that, similar to previous works^, 
both T and T e s act as a mass for the quantum fluctuations. At T = and T e g = 0, x~k vanishes at low energies as ~ |Q|. 
However, for T + or T^r + 0, Xk ^ as Q — > 0, and the theory effectively becomes classical. 



2. RG analysis in eigenmode basis 

We now perform a renormalization group analysis of the effective action, Eq. dDlOb . Following a similar analysis as in the 
main text, we find the following set of RG equations to lowest (linear) order in the quartic coupling, 



dlnb 

d\'d{b) 
dlnb 
dT(b) 

dlnb 

duj(b) 

dlnb 
dA(b) 
dlnb 



3T eff (b), (D15) 

2v d (b), (D16) 

-r(b), (D17) 

-u t (b)+0([Uj] 2 ), (D18) 

2A(b) + 3umf{T eS {b), rib)). (D19) 



The function f(T e ff(b), r(b)) is given by 



2tt 

where Xk(p) = ~X'r(p)/\X'r'(p)\ 2 - Solving these equations, we arrive at 



f(T eB (b),T(b)) = iK 2 J ^XKt 



T ee (b) = T eB b\ (D21) 

v d (b) = v d b\ (D22) 

Tib) = Tb~\ (D23) 

Uiib) = Uib-\ (D24) 

Jr»ln b 

dxe- 3x f(T eS e x \ re - *)]. (D25) 
o 



Tib) = Tb\ (D26) 



An analogous renormalization group analysis in the equilibrium case gives 

3 

T(b) = Tb-\ (D27) 
Uiib) = Uib-\ (D28) 

Jr»ln b 
dxe~ 3x fiTe 3x , re - *)] . (D29) 
o 
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We now recompute the correlation length in the eigenbasis. The scaling stops when A(bi) ~ 1. In the quantum disordered 
regime, T and T^g is small enough so that T{b\) <*c 1 and T e g{bx) <sc 1. In this case, the correlation length of the system can be 
obtained by setting T — and T e g = 0. One then obtains 

^ 2 = A + 3mi I dxe~ 3x f(0, re~ x ) =: r. (D30) 
Jo 

The condition for the occurrence of the regime is obtained by requiring T(b\) < 1 where b\ = r Then the condition reads 

T < r 3 ' 2 , T eB < r 3 ' 2 . (D31) 

The non-universal shift to the bare mass in Eq.( lD30l l indicates that we are above the upper critical dimension and that u\ is a 
dangerously irrelevant operator. 



When the inequality in Eq.( ID3 It is violated, the critical fluctuations become classical in nature and the system enters the 
quantum critical regime. Here we compute the leading corrections to the correlation length in this regime. We scale up to 
either T(b*) = 1 or T e ff{b*) = 1. Recall that the solution for A(b) in equilibrium is given by Eq.( |D29t . The correction to the 
zero-temperature distance from critic ality, r, is then given by 



/-»ln b* 

5r = 3«i dxe- 3x [f(Te 3 \ re - *) - /(0, re - *)]- 

Jo 



(D32) 



If we define Af :- f(Te , re x ) - f(0, re *), we have 



Af = iK 2 



I 



dQ. 
2^ 



(D33) 



and 



a^(i,q) = -in 



coth — 
2T 



sgn(O) 



t+ 1 



(A+ l) 2 + (Qr + Q) 2 /4 



(D34) 



The Q-integral gets most of the contribution from -2T < Q. < 2T. So we can cutoff the integral there, and use coth(f2/2r) 
277 £1 Then we obtain 



4*rr(i,n)*-i(2r-|n|) 



1 +T 



1 + (Qt + Q) 2 /4' 



(D35) 



where we dropped A from the denominator since it remains small during the flow. The Q-integral can now be done trivially. 
Expanding the result for small T, we get 



a/= l(i +T )r 2 +0(r 3 ). 



(D36) 



To evaluate the integral up to b*, we first do a change of variables, x — > T(b) -: y, which makes the measure dx = dy/3y. The 
upper limit of the integral is 1, and we extend the lower limit down to 0. Then the correction is given by, 



8r 



Jo 3y y 



T 



l/3\ 



(D37) 



We then obtain 



8r = -4;T 1 + — T 



3t 



'1/3 



(D38) 



We now move on to the nonequilibrium case. Recall that the solution A(b) for the nonequilibrium case is given by Eq.( ID25l ). 
The correction to the zero-temperature distance from criticality, r, is then given by 

dxe- 3x [f{T e $e 3 \ re~ x ) - f(0, Te~ x )]. (D39) 

o 
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We then obtain 



A^(l,Q) = -/|Q|r 



1 +(Qt + Q) 2 /4' 



(D40) 



where we have dropped the drift and the mass terms from the denominator since they remain small during the flow 8 . Once again, 
the f2-integral is appreciable only for |£2| < 2T e g. In this region, we have 



j(]js\ „ ?eff 2 

\\Q,\ )* m| n 



Performing the f2-integral up to 2T^, we then obtain to lowest order in T e g, 

A /-4< ff+ o(r e 3 ff ). 



(D41) 



(D42) 



The mass correction is then given by 



Performing the integral, we obtain 



f l dyT e *J (T e A m ) 

6r = 3 Ul - A/ y, T 

Jo 3y y \ \ y I ) 



5r = ^i Tr 4/3 
0r n 3 ' eff ■ 



To summarize the results, the correlation length in equilibrium reads 



9 ui^l, 3t 



'1/3 



while out of equilibrium, one obtains 



r 1 = r + — - T 

?neq ^3 eff 1 



3U]T 



l T T 4/3 



(D43) 



(D44) 



(D45) 



(D46) 



This is in exact agreement with the results obtained in the main text (c.f. Eqs. (|90]l,(l92l)). 
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